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ABSTRACT 

The mass symmetry between the two species in electron-positron (pair) plasmas has inter- 
esting consequences for collisionless magnetic reconnection because the Hall term, which plays 
a crucial role in supporting fast reconnection in electron-proton plasmas, vanishes. We perform 
kinetic simulations of pair reconnection in systems of various sizes, show that it remains fast, 
and identify the reason why this occurs. For sufficiently large systems a Weibel-like temperature 
anisotropy instability develops in the outflow from the X-point that causes the current layer to 
broaden and form a Petschek-like open outflow. We discuss the parameter regimes in which pair 
reconnection should be fast and the implications for astrophysical pair plasmas. 

Subject headings: magnetic fields — methods: numerical — plasmas 



1. Introduction 

The energy density of magnetic fields can be 
significant, and even dominant, in certain as- 
trophysical systems and so the question as to 
which processes act as sinks and convert that 
energy into other forms naturally arises. Obser- 
vations (including in situ measurements) show 
that magnetic reconnection, in which topological 
changes in an embedded field lead to increases 
in a plasma's kinetic and thermal energy, plays 
such a role in the terrestrial magnet osphere, the 
solar wind, and the so l ar corona. §Mozer et 



20021 : \Phan et al\ 120061 : \Masuda et a/.lll994f ) The 
natural extension to other astrophysical locations 
has led to the consideration of the role of recon- 



nection in accretion disks bEardlev and Liahtman 



19751; \Miller and S tone 1997), gamma-ray bursts 
( Drenkh ahn and Svruitl 120021) and pulsar winds 
( Cor oniti\\199(x I MichelW 1994 : 1 Zenitani and Hoshino 
20011 ). Unlike the heliosphere, positrons, rather 
than protons, can comprise the positively charged 
species in the plasmas of these systems. The equal- 
ity of the electron and positron masses introduces 



a symmetry that can have a potentially large effect 
on the development and role of reconnection. 
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Parker and Sweet \Parker\ Il957c \Sweet\ 11958) 
proposed the first model of electron-proton recon- 
nection, a magnetohydrodynamic (MHD) formula- 
tion that attempted to explain the fast timescales 
associated with solar flares. However contempora- 
neous observations suggested that magnetic flux 
reconnected much faster than the model predicted 
(on timescales of minutes, not weeks). Relaxing 
the assumptions of classical resistive MHD by in- 
voking an anomalously large and localized resis- 
tivity, perhaps driven by turbulence, resolved the 
theoretical discrepancy, but no model of the pro- 
cess garnered widespread support. 

Later work focused on Hall MHD, an exten- 
sion of MHD that includes the decoupling of the 
motions of protons and electrons at small spa- 
tial scales — on the order of the ion skin depth, 
di = c/uu p i, where uo pi is the proton plasma fre- 
quency. The GEM Challenge, a collaboration 
within the reconnection simulation community, 
showed that various algorithms including the Hall 
term (the simplest new term arising in Hall MHD) 
all reconnected at approximately the same rate: 
V[ n « O.Iva, where ^i n is the speed at which 



1 



plasma flows into the X-point and va is the lo- 
cal Alfven speed iBirn et all 12001 ). A resistive 
MHD code simulating the same problem recon- 
nected flux at a rate several orders of magnitude 
slower. 

Although the GEM Challenge demonstrated 
that inclusion of the Hall term is a sufficient re- 
quirement for fast reconnection, the question of 
whether it is a necessary requirement remains 
open. Pair reconnection provides an excellent test 
because one effect of the mass symmetry is the 
disappearance of the Hall term from the govern- 
ing equations. If it is necessary for fast reconnec- 
tion then pair plasmas will reconnect slowly, on 
Sweet-Parker timescales, and magnetic reconnec- 
tion will be an inefficient mechanism for releasing 
the magnetic energy in pair plasmas. 

Despite this interesting limit the simulation 
community has only r ecently begun to consider 
pair r econne ction. \Bessho and Bhattachariee 
(|2Q05l . l2QQ7h observed fast pair reconnection 
in a small system and attributed it to a lo- 
calized resistivity-like effect caused by the off- 
diagonal components of t he p ressure tensor. 
Dauahton and Karimabadn ([20071 ) did a careful 
examination of pair reconnection in large systems 
and arrived at similar conclusions for the early, 
impulsive phase of reconnection. Over longer 
periods they suggested that generation of large 
numbers of magnetic islands served to localize 
the current layer a nd support fa st reconnection. 
Zenitani and Hoshind (l200lh , \jaroschek et al. 
( 2004 ). and \Fuiimotd (|2006l ) explored other as- 
pects of pair reconnection but, although they 
observed fast reconnection, did not advance any 
causal explanations. 

We also find that fast reconnection occurs in 
pair plasm as. However, we reach differ e nt con clu- 
sions from \Dauahton and KarimabaZi ( 20071 ) on 
the role of magnetic islands. In the initial phase 
of electron-positron reconnection many magnetic 
islands typically form, especially if the ratio of the 
initial width of the current layer to the length of 
the system is very small. These islands have a 
transient influence on the development of recon- 
nection until they convect downstream. Secondary 
islands also form later in time, and in greater num- 
bers than in electron-proton reconnection. But, 
we find that these islands do not seem to control 
the overall length of the current layer and there- 



fore the rate of reconnection. Based on our simu- 
lations we instead propose that the outflow jet in 
large pair plasma systems opens as a result of tur- 
bulence. We show that as the system size grows 
the current layer lengthens, just as Hall- less Sweet- 
Parker theory predicts. Interestingly the recon- 
nection rate remains constant during this change 
for small systems because other parameters, pri- 
marily the speed of the outflow jet, also vary. 
This cannot continue for arbitrarily large boxes 
because of the Alfvenic limitation on the outflow 
speed. Instead, once the system is big enough 
(and the current layer long enough) an electromag- 
netic pressure-anisotropy instability similar to the 
Weibel mode begins to grow in the outflow jet. 
The outflow jet becomes strongly turbulent and 
the outflow layer b roadens downs tream as pro- 
posed by Petschek ^Petschekl fl964). We present 
a simple analytic argument that suggests that as 
a result of this turbulence the rate of reconnection 
should take on a constant value, independent of 
system size. 

In Section[2]of this paper we describe the details 
of our computations. Section [3] outlines the role 
of the Hall term and how it disappears in pair 
plasmas. We discuss the simulations in Section [4] 
and the instability and its consequences in Section 
[5l Finally, Section [6] explores some implications of 
our findings. 

2. Computational Details 

We use p 3d, a massively p arallel particle-in-cell 
(PIC) code \Zeiler et aZ.ll2QQ2h to perform our sim- 
ulations. As in other PIC codes we divide the 
computational domain into cells, track the electro- 
magnetic fields only on the gridpoints, and allow 
the particles to move freely. During a timestep 
the Boris algorithm advances the Lorentz equa- 
tion of motion for each particle: the electric field 
E provides an acceleration for half a timestep, 
the magnetic field B rotates the velocity vector, 
and E accelerates for the second half- timestep. 
An explicit trapezoidal-leapfrog method employ- 
ing second-order spatial derivatives in Faraday's 
and Ampere's Laws advances the fields in time. 
We use uniform, square computational cells that 
force V'B = 0. However our discretization scheme 
leads to a discrepancy between V • E and 47rp. To 
address this problem we use a multigrid algorithm 
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to solve Poisson's equation for a correction term 
to E. Our particles do not directly interact and, 
in particular, no electron-positron annihilation oc- 
curs. 

To help elucidate the underlying physics we 
write the code equations in normalized units. 
Masses are normalized to the ion mass m^, mag- 
netic fields to the asymptotic value of the reversed 
field, and the density to the value at the center 
of the current sheet. Other normalizations de- 
rive from these: velocities to the Alfven speed va, 
lengths to the ion inertial length c/u P i, times to 
the inverse ion cyclotron frequency , and tem- 
peratures to rriiV 2 A . In electron-positron plasmas 
m e /rrii = 1. 

In our coordinate system the inflow and outflow 
for an X-point lie parallel to y and x, respectively. 
The reconnection electric field is parallel to z. In 
the simulations presented here we assume out-of- 
plane derivatives vanish, i.e., d/dz = 0; this choice 
eliminates any structure in the z direction. The 
reasonableness of simulating pair plasmas with 
this restriction remains unknown. For electron- 
proton plasmas, a comparison of 2D and 3D recon- 
nection simulations suggests that the qualitativ e 



features remain unchanged (Hess e et ali I20011 ). 



However for many instabilities that can have large 
effects on the out-of-plane structure, e.g., the drift- 
kink, the onset threshold varies with the mass ra- 
tio, and so further simulations are needed to ex- 
plore their importance in the m e /mi = 1 limit. 

The initial e quilibrium c onsists of two Harris 
current sheets ^Harris} 1 19621 ) superimposed on a 
ambient population of uniform density. The re- 
connecting magnetic field is given by B x /Bq = 
tanh[(j/ - L y /4)/wo] - tanh[(j/ - 3Ly/4)/wo] - 1, 
where wo and L y denote the half-width of the 
current sheets and the box size in the y direc- 
tion respectively. This configuration allows us 
to use fully periodic boundary conditions. Both 
species have the same uniform initial tempera- 
ture, T e = Ti = 0.25. Pressure balance uniquely 
determines the density profile, aside from a uni- 
form background population that can have arbi- 
trary density (here = 0.2); in this equilib- 
rium n(y Ly/4) = n(y = 3L y /4) = 1.2. At 
t — we perturb the magnetic field — B x /Bq 
varies with system size but is always < 0.07 - 
to seed X-points at (x, y) = (L^/4, 3L y /4) and 
(3^/4, L y /A). 



We set the speed of light (in normalized units) 
to 5. The spatial resolution is such that there are 

> 4 gridpoints per inertial length and « 1 per De- 
bye length. The Courant condition determines the 
particle timestep; we substep the advancement of 
the electromagnetic fields. A typical cell contains 
~200 particles and our largest simulation follows 

> 10 9 particles. All of our simulations conserve 
energy to better than 1 part in 200. 

3. The Hall Term 

The current theory of electron-proton reconnec- 
tion emphasizes the role of the Hall term in the 
generalized Ohm's Law. By combining the col- 
lisionless fluid momentum equations for electrons 
and a singly charged positive ion species of mass 
rrii we can write 
1 



(l+/i)E: 



/i 



v x B 



1 



/i 



nec 
ne 



J x B 



dJ 



Jv + vJ 



1 l-/i 

ne 1 + \i 
(1) 
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where a = m e /mi, v = (m e v e + miVi)/(m e + mi) 
is the velocity in the center of mass frame, P is the 
pressure tensor, and we assume quasi- neutrality, 
rii = n e = n. For electron-proton plasmas a w 0. 

By scaling equation [1] one can show that 
the first term on the right-hand side dominates 
at large lengthscales. The MHD model, and 
hence Sweet-Parker theory, discards the remaining 
terms. The second term on the right-hand side is 
the Hall term; it becomes important for length- 
scales equal to or smaller than With the addi- 
tion of the Hall term the principal wave mode on 
these scales changes from the Alfven wave of MHD 
{uj ex k) to the dispersive (uo (x k 2 ) whistler wave. 
The localization of the current layer along the out- 
flow direction due to the action of dispersive waves 
has been proposed as the mechanism that facili- 
tates fast reconnection in electron-pr oton plasmas 



\Birn et aZ.ll200ll : \R,oaers et a/.ll200lf ). The terms 



proportional to the diagonal portion of the pres- 
sure tensor can be important for systems with 
large guide fields (since they are the source of ki- 
netic Alfven waves), while the off-diagonal terms 
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and final inertial term matter only on electron 
scales, d e = c/ou pe , and likely play only a small 
role outside the immediate vicinity of the X-point 
and magnetic separatrices. 

In pair plasmas \i = 1 and equation [T] simplifies 
to the form derived in \Bessho and Bhattacharjee 



(2005): 



E 



- v x B 

c 

^~ V • (Pe 

2ne 

m e |~<9J 
~gi^ 



-Pi) 

V • (Jv + vJ) 



(2) 



The mass symmetry eliminates the Hall term, and 
with it the usual dispersive modes, including the 
whistler. If we make the reasonable assumption 
that v vanishes at a steady-state X-point, equation 
[2] becomes 



E = 



(3) 



Hence a non-zero reconnection electric field (E z ) 
requires an asymmetry between the two species, 
P e 7^ P^. In our quasi- two-dimensional system the 
balance must come from the off-diagonal compo- 
nents of the pressure tensors, which describe the 
transport of momentum away from the X-point. 

4. Simulations 
4.1. Overview 

Figure [1] shows the steady-state behavior of a 
typical simulation of pair reconnection in a large 
domain. The top panel displays magnetic field 
lines superimposed on the component of the elec- 
tron velocity in the x direction; the positron ve- 
locity (not shown) is basically identical. The solid 
and dashed portions of the field lines indicate pos- 
itive and negative signs of B x , the reconnecting 
component of the field. The tension in the highly 
bent field lines to the left and right of the X-point 
(x = 0, y = 0) accelerates the plasma downstream 
until it reaches the local Alfven velocity (~ 1 in 
our normalization). A more gradual inflow with 
speed ^0.1 (not shown) replenishes the plasma. 

Panel (b) shows the out-of-plane component of 
the magnetic field B z . The quadrupolar B z pat- 
tern that is one of the signatures of Hall reconnec- 
tion is not present. Instead B z becomes large and 



highly structured farther downstream from the X- 
point where the current layer opens up. We dis- 
cuss the turbulence, the instability responsible for 
it, and its linkage to the structure of the outflow 
jet later. 

Although E y has a structure similar to B z in 
the turbulent outflow, along the thin central por- 
tion of the current layer it exhibits a narrow bipo- 
lar feature. Panel (c) shows a vertical cut of E y 
averaged over the region —50 < x < 50. The cen- 
tral structure has a width of « 2d e and breaks the 
symmetry between the electrons and positrons. It 
grows from the ambient noise in all of our simula- 
tions, regardless of box size, but the sign (whether 
electrons or positrons accelerate toward the cur- 
rent layer) appears to arise spontaneously. In 
this case electrons are accelerated toward, and 
positrons repelled from, the center, which leads 
to an asymmetry in the respective particle dis- 
tribution functions. The bipolar signal appears 
much earlier in the simulations than the temper- 
ature anisotropy instability we discuss later, and 
it appears in runs for which this other instability 
is absent. Instead, the counter-streaming inflow- 
ing plasma likely triggers a two-stream instability 
at the current layer in which the polarity of E y is 
determined randomly. 

To maintain this field the species' number 
densities exhibit a slight asymmetry (i.e., quasi- 
neutrality is slightly violated). Panel (d) displays 
the distribution of v y for the electrons (green) and 
positrons (blue) within the box \x\ < 50, \y\ < 2 
and shows the effects of this field. The v z distri- 
butions in panel (e) are mirror-symmetric due to 
the effect of E z (the reconnection electric field) 
on the oppositely charged electrons and positrons. 
The v x distributions are the same for both species 
because E x is small. They are plotted in black. 

In some of our preliminary simulations the long 
initial current layers produced a larger number of 
islands that played a transient role before being 
swept away. We find that varying the thickness 
of the initial layer affects the initial rate of island 
production, with a broader layer producing fewer 
islands. At the time of Figure [T] these initial is- 
lands have all disappeared downstream. However, 
we also find that secondary islands self-generated 
on the current layer remain a common occurrence 
in pair reconnection and are much mo re prevalent 
than in reconnection with m e <C mi (Sha y et al 
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20071 ). Two such islands (at x « 20 and x « -80) 
appear in Figure [H These remnant islands, how- 
ever, are not responsible for the turbulent fluctu- 
ations in B z in Figure [H The island at x ~ 20, 
occurs in a region with no turbulence and has no 
£^ signature while the other, at x ~ —80, lies in 
the turbulent outflow where any independent B z 
signature (if it exists) is difficult to distinguish. 

4.2. Scaling of the Reconnection Rate 

Fast reconnection ultimate ly arises from the lo- 
caliza tion of the current layer \Biskamp and Schwarz 
l200ll ). This can be seen from a straightforward 
analysis of the fluid continuity equation 

^+V.(nv)=0 (4) 

In a steady state the time derivative vanishes, leav- 
ing a term that the divergence theorem transforms 
into a surface integral. For one quadrant of a rect- 
angular domain with length (parallel to x) 2A and 
width (parallel to y) 25 centered on the X-point, 
the surface integrals imply that 



^in = ^out 



™out 



A J \ m n 



(5) 



where the subscripts "in" and "out" refer to 
inflowing and outflowing plasmas, respectively. 
Both the outflow speed, which approaches the lo- 
cal Alfven speed, and the density ratio n ou t/^in, 
which does not radically change from its initial 
value, are 0(1). If, as seems a reasonable assump- 
tion, the inflow scale 5 equals an intrinsic kinetic 
scale-length such as d e and remains roughly con- 
stant as the system size changes then v in scales 
inversely with the outflow dimension A. In par- 
ticular, if, as in the Sweet-Parker picture, A scales 
with the size of the system, V[ n should decrease for 
larger and larger simulation domains. In electron- 
proton reconnection the Hall term keeps A fixed, 
regardless of the system size, and maintains a 
fast reconnection rate Vi n /v out ~ 5/ A ~ O(0. 1) 
\Shav et a/.lll999i l2007h . 

If the Hall term is the only mechanism that 
can localize the current layer then pair reconnec- 
tion should feature current sheets that scale with 
the system size. In order to test this behavior 
we began with a relatively small system and then 
doubled the size repeatedly. Table [1] contains the 



parameters of the four runs we discuss in this pa- 
per; Figure [1] comes from run d. Because of com- 
putational constraints we doubled L x but not L y 
between c and d, and so the latter's aspect ratio 
differs from the other three runs. Earlier work 
showed that if the system is large enough to reach 
a steady state the reconnection rate in a double 
tearing mod e configuration is independent of the 
aspect ratio \Shav et aZ.lll999h . 

Figure [2] shows the steady state current layers 
for the four runs. In the smallest simulations (pan- 
els (a) and (b)) the current layers lengths approach 
the system size. But, when the system size dou- 
bles between panels (b) and (c) the current layer 
does not; it lengthens moderately, but the turbu- 
lence evident at both ends of the layer broadens 
it along the y direction and prevents the narrow 
portion of the layer from expanding further. A fi- 
nal lengthwise doubling of the system, panel (d), 
barely increases the length of the current layer. 

In Figure [3] we show the reconnection rates 
as a function of time. To measure this rate we 
first integrate the magnetic flux between the ini- 
tial X-point and the large magnetic island formed 
by the reconnected field. (The large island is a 
consequence of our periodic boundaries and is not 
shown in the figures.) We check for the growth of 
other X-points, but find that the initially seeded 
one remains dominant throughout the runs, de- 
spite the formation of secondary magnetic islands. 
Figure [3] shows the (unsmoothed) temporal deriva- 
tive of the reconnected flux. The vertical lines 
mark the times of the snapshots of v ez displayed 
in Figure [21 All of the runs exhibit essentially the 
same behavior: after an initial quiescent period 
the reconnection rate rises to an asymptotic value 
of ~ 0.1 and remains there. Given enough time 
the separatrices bounding the magnetic island on 
the simulation's other current layer approach the 
X-points shown in Figure [2] and we halt the run. 
The first effects of this interaction can be seen in 
the slight increases in the reconnection rates no- 
ticeable at the end of each panel in Figure 03 

So, despite the lack of a Hall term, pair re- 
connection remains fast when the system size in- 
creases. A more quantitative description of the re- 
connection can be found in Table [21 which lists the 
various quantities from equation [5] for each run. 
The 

^in.sim column contains the inflow speed mea- 
sured in the simulation while V[ n ca i c is the value 
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calculated from the other data using equation [5j 

For our smallest domain, run a, 2A w 35 ap- 
proaches L x /2 = 50. (The periodic boundaries 
in our simulations restrict the maximum current 
sheet length to L x /2. In practice A never reaches 
this maximum because the transition from out- 
flow to island has a finite scale.) Even though 
the current layer stretches to the system size in 
a Sweet-Parker-like manner, reconnection remains 
fast (vin ~ O(0. 1)) because the box is relatively 
small. In fact the outflow velocity does not even 
reach the Alfven speed. Thus, exploring the scal- 
ing of reconnection in pair plasmas requires much 
larger simulation domains than in normal recon- 
nection. In a simulation domain of this size (mea- 
sured in ion inertial lengths) the ion outflow ve- 
locity would have reached the full Alfven speed. 

After doubling the box size (run b) , the current 
layer length A also increases by a factor of two. 
However the larger box contains more room for 
the acceleration of the outflow and so v out also in- 
creases. The densities, both upstream and down- 
stream, as well as the thickness of the current layer 
remain relatively constant between these two runs, 
and all of our other runs. The net effect of the 
variation of A and v out is to leave the reconnec- 
tion rate unchanged. Since the outflow speed has 
an upper limit, the Alfven speed, this progression 
cannot continue indefinitely. 

After the next doubling (run c) the outflow ve- 
locity reaches its expected Alfvenic limit. Since 
the increase was less than a factor of 2 (1.3/0.8 ~ 
1.6) a Sweet-Parker-like doubling of A would be 
accompanied by a decrease in the reconnection 
rate. Instead of doubling, however, A only in- 
creases by a factor of 120/80 « 1.5 and the recon- 
nection rate stays approximately constant. Since 
the Hall term cannot halt this increase, some other 
process must be at work. 

As a check we again doubled the length of 
the box (but not its width) for run d. Unlike 
in the previous doublings the outflow speed can- 
not increase further. Its maximum value is the 
Alfven speed based on the downstream density 
and the upstream magnetic field; in these simu- 
lations max(^out) ~ 1.3. Absent any other localiz- 
ing effect the current layer length should approx- 
imately double and the reconnection rate should 
halve. Instead the current layer length and the 
reconnection rate remain basically unchanged. In 



Section [5] we argue that the length of the current 
sheet is limited by the development of an instabil- 
ity driven by a temperature anisotropy. The in- 
stability develops in the outflow from the X-point 
before eventually growing strong enough to open 
the layer. 

5. Anisotropy Instability 

In sufficiently large simulations the outflow 
from the X-point is subject to a temperature 
anisotropy instability. The top panel of Figure 
2] shows a view B z from run c at t = 450, soon 
after the instability develops. In the second panel 
we plot two of the diagonal components of the 
positron temperature tensor, T xx and T yy , along 
the line y = — 1. (The influence of the E y feature 
shown in Figure [T] makes the anisotropy somewhat 
weaker exactly on the symmetry axis.) At the X- 
point (x = 0) the plasma temperature is isotropic 
and equal to the initial value of 0.25. As the 
plasma moves away from the X-point T xx increases 
sharply while T yy remains nearly unchanged, lead- 
ing to an anisotropy that reaches a peak value of 
«4. At \x\ « 50 the disturbance in B z becomes 
visible, indicating onset of an instability. Further 
downstream the B z pattern strengthens and the 
anisotropy simultaneously weakens. Movies of B z 
(not shown) demonstrate that the disturbances 
propagate with the mean flow of the plasma in 
the outflow jet. 

In order to investigate the source of the 
anisotropy we show the 1-D velocity distributions 
taken between the dashed lines in panel (b) and 
— 2 < y < — lin the final three panels. Two pop- 
ulations intermingle in this region, the inflowing 
plasma with v x ~ and a jet with v x < moving 
outward from the X-point. As can be seen in panel 
(c) the superposition of these populations results 
in a wide, hot distribution of v x . Meanwhile, 
since the sample box is upstream of the bipolar 
E y signal, the v y distributions remain essentially 
unchanged from their initial Maxwellians. A mix- 
ture of the inflowing ambient plasma and plasma 
accelerated by the reconnection electric field form 
the v z distributions. The final result is a plasma 
with T x > T y « T z . 

Plasmas with hot bimaxwellian distribution 
functions are typically stable to electrostatic per- 
turbations (the thermal spread stabilizes the elec- 
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trostatic ion- ion streaming instability), but unsta- 
ble to electromagnetic instabilities. One possibil- 
ity is that the magnetic perturbations in Figures 
[D and [H result from the firehose instability. If that 
were the case the perturbed motion of the plasma 
in the z direction with k = fc^x would produce per- 
turbations in B z . We can rule out this instability, 
however, since along the symmetry axis (y = 0) 
motion of the plasma in z would not produce a 
perturbation of B z since the equilibrium field B x 
is zero in this region. The simulations clearly 
have large perturbations of B z centered at y = 0. 
Instead the anisotropy is sufficiently large that 
the magnetic field in the region around y = is 
negligible and the instability is a Weibel-like mode 
that self generates the magnetic field perturbation 
by separating the electrons and positrons moving 
in opposite direction to create a current in the 
x-direction that varies with y. To demonstrate 
the viability of the instability we begin with a tri- 
maxwellian plasma described by the distribution 
function 



the growth rate of 



n 



\2TTk B T a 



1/2 



exp 



m 



2nk B T a 



(6) 

with, in general, T x ^ T y =^ T Zl and look for insta- 
bilities with B z and E x vary ing as exjpjik^y — tout) . 
Following the method of \Krall and Trivelpiece 
(1986), the instability condition for a pair plasma 
is 



1 > 



Setting the plasma dielectric function to zero, we 
obtain the dispersion relation relating the growth 
rate 7 = —iuj and the wavenumber: 

V 2 +7 2 +2^(l-^)=2^| CZ (C) (8) 



y / 

where the thermal speed v t h enters through 



k y yj2T y /m ky v t h 
and Z(Q is the plasma dispersion function 



(9) 



Z (o4f^* do) 



J -00 t - c 

In the limit \(\ ^> 1 equation [8] has a solution for 



7 ~ k y vth,x 



k 2 y c 2 + 2uo 2 p 



(11) 



where ^th,x = \/^T x /m. Thus, the growth rate 
asymptotes to 7 ~ v t h,x/de for k y d e > \/2. Taking 
ICI ^> 1 is only possible if the anisotropy is strong, 
T x /T y ^> 1. Equation [TT] can also be derived in 
the fluid limit by considering the stability of a cold 
plasma beam. 

Strictly speaking these results apply only to the 
case of a homogeneous infinite plasma, while in 
the case under consideration here the instability 
develops within a thin (width 26 w 5d e ) current 
layer. To further explore the behavior of this in- 
stability, particularly the behavior of its growth 
rate, we again turn to simulations. First, in or- 
der to test our code's ability to reproduce the re- 
sults of equation [51 we consider a homogeneous 
electron-positron system with a bimaxwellian dis- 
tribution (for both species), no initial magnetic 
field, and a constant particle density, n = 1.2. We 
seed the system with the theoretical fastest grow- 
ing wavenumber (with k x = 0) and monitor the 
growth of the instability. In panel (c) of Figure [5] 
we plot the growth rates calculated from the simu- 
lations (stars) compared with the values predicted 
by equation [8] (solid curve). The two are in near 
perfect agreement. 

In order to study the effect of the narrow cur- 
rent layer on the development of the instability we 
next turn to a more complicated system. We begin 
with a Harris sheet equilibrium similar to that dis- 
cussed in Section [2] except that one component of 
the temperature, T X: increases along with the den- 
sity in the center of the current layer. Note that 
T x does not affect force balance in the y direction 
and therefore does not affect the initial equilib- 
rium. For our usual box sizes this equilibrium is 
also unstable to the tearing mode, so we consider 
narrow domains (e.g., L x = 4, L y = 50) where 
the tearing mode is stabilized but the temperature 
anisotropy mode is not. In panel (d) of Figure [5] 
we show the growth rates (marked by the stars) 
for the instability in this geometry (again with 
k x = 0) compared with the homogeneous growth 
rate (solid curve). In general the growth rates are 
only reduced around 20 % from the homogeneous 
values. Thus, neither the finite geometry of the 
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current layer nor the finite magnetic field in the 
layer significantly affect the growth of the insta- 
bility, confirming that the Weibel instability can 
exist in a magnetically supported current layer. 

In panels (a) and (b) of Figure [5] we show the 
nonlinear development of the instability in a 2-D 
homogeneous system with no initial magnetic field 
and T x /T y = 4. The system is seeded with a mag- 
netic perturbation corresponding to the fastest 
growing mode that dominates the structure at 
early time. At late time the dominant mode de- 
velops a finite value of k x (k x = 0.25, k y = 1.0), as 
also occurs in the reconnection simulations. The 
amplitude of the late-time perturbations of B z is 
about 0.5; compared with the turbulence seen in 
the reconnection simulations these perturbations 
are at somewhat longer wavelength and have mod- 
estly larger amplitudes. These differences may re- 
flect constraints on the turbulence due to the fi- 
nite width (in the y direction) of the temperature 
anisotropy in the reconnection simulations. 

For the full reconnection simulation a typical 
temperature anisotropy in the outflow is w 3 (see 
Figure [4j) , which corresponds to a typical growth 
rate of 7 ~ 0.2. (The density in the reconnection 
outflow is w 4 times smaller than in the insta- 
bility simulations, which, since 7 oc uo pi implies a 
factor of 2 decrease in the growth rate.) In plasma 
traveling at the typical outflow velocity of v out « 1 
the lengthscale associated with one e-folding of the 
instability is t ~ v ut/7 = 5. A current layer of 
half-length A = 60 can then fit A/£ « 12 growth 
periods within it, more than enough to amplify 
any initial fluctuations to non-linear amplitudes. 

A fundamental question is whether the dynam- 
ics of the Weibel instability scale in a manner such 
that the resultant length of the current layer can 
produce reconnection rates insensitive to the size 
of the computational domain. A complete model 
should incorporate the development of the temper- 
ature anisotropy in the outflow from the X-point 
by using the full dispersion relation of equation [HI 
We consider instead a cruder model in which we 
ignore the anisotropy threshold. That is, we take 
the plasma near the X-point to be cold so that the 
Weibel mode is unstable as soon as the anisotropy 
develops downstream. Strictly speaking this as- 
sumption does not apply near the X-point and 
will lead to an overstatement of the mode's growth 
rate. But since the principal effects of the instabil- 



ity manifest themselves in the downstream region, 
the key result derived from this simpler approach 
— the independence of the reconnection rate on 
the system size — should remain valid in a more 
complete theory. 

In its most primitive form, neglecting onset 
thresholds, the growth rate of the Weibel is given 
by 7 = Vth,x/de- I n Figured! we show that the ef- 
fective electron thermal speed in the x-direction 
is related to the relative drift between the ba- 
sically stationary upstream plasma crossing into 
the outflow jet and the high speed particles ac- 
celerated from the X-point. Hence Vth, x is sim- 
ply one-half the drift speed of those fast particles, 
which implies 7 = v x /2d e . As discussed previ- 
ously, the amplitude of the turbulent fluctuations 
is controlled by the convective amplification of the 
initial electromagnetic perturbations. The num- 
ber of e-foldings between the X-point and an po- 
sition x downstream is therefore given by 

JO v x 

But the growth rate 7 is itself proportional to v x , 
so v x drops out of the calculation and 

T = x/2d e . (13) 

If we assume that To is the number of e-foldings 
required to amplify initial fluctuations to the val- 
ues of order unity required to broaden the current 
layer, we obtain the length of the current layer 
A = 2Tod e . The resulting rate of reconnection 
can then be calculated from equation [5] (ignoring 
density changes and any variations in the current 
layer width) to be 

^in ~ v A /T . (14) 

The rate of reconnection resulting from the broad- 
ening of the current layer by Weibel generated tur- 
bulence should therefore remain fast, independent 
of system size, and independent of kinetic scales, 
as we see in our simulations 

6. Discussion 

We have shown that reconnection in a pair 
plasma remains fast despite the absence of the 
Hall term. For small systems, up to 100<i e along 
the outflow direction, the length of the current 
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layer that develops during reconnection scales with 
the size of the system, as expected based on the 
classical Sweet-Parker model with no dispersive 
waves. When the domain size exceeds ~ 200d e 
a Weibel-like instability develops that forces the 
development of a turbulent outflow jet and lim- 
its the length of the narrow current layer, allow- 
ing reconnection to remain fast. Note that, de- 
spite being big enough to show the transition to 
system-size independent reconnection, our simula- 
tion boxes remain small in astrophysical terms. A 
scale of L = 10 3 inertial lengths, modestly larger 
than our biggest boxes, corresponds to 



L ~ 5 x 10 s 



1 cm 



1/2 



km 



(15) 



The length A w 70 d e and width S ~ 5 d e of the 
current layer are even smaller still. The wide sep- 
aration of scale between the system size and the 
scale where the important reconnection physics 
occurs is a well-known phenomenon in magneto- 
spheric physics where in situ satellite measure- 
ments have studied the inner scales near the X- 
point \Mozer et aZ.ll2QQ2h . 

Although fast pair reconnection has been noted 
in previous work, the underlying cause was at- 
tri buted to different mechan i sms. In particu- 
lar \D<mghtm} L j^^ (|2QQ7I ) suggest that 
the generation and ejection of magnetic islands 
play an important role by continually disrupting 
the growth of a current layer that would other- 
wise scale with the system size. Magnetic islands 
form in our simulations, but appear to play little 
or no role in maintaining the length of the cur- 
rent layer. Figure [6] shows snapshots at intervals 
of fldSt = 25 of the out-of-plane electron velocity 
for run <i, beginning at Q C it = 600. During this 
period the reconnection rate remains steady and 
fast. In the first four panels the instability de- 
scribed in section [5] does not have enough room to 
grow and the current layer lengthens, in a process 
analogous to what occurs in our smaller simula- 
tion domains. In the second column the current 
layer grows long enough for the Weibel-like insta- 
bility to saturate and stop further lengthening of 
the narrow portion of the current layer. Although 
several islands form and travel downstream during 
this time, they do not appear to have any signif- 
icant effect on the length of the current layer or 
on the system's reconnection rate (see Figure [3]). 



However, several factors may play a role in island 
generation in numerical simulations including the 
initial width of the current sheet, the strength of 
the initial perturbation, the number of simulated 
particles, and the size of the grid. The relative 
importance of these various effects, and what they 
imply about island generation in real plasmas, re- 
mains unknown. 

A surprise of the simulations of pair reconnec- 
tion shown here and in earlier work is that, just 
as in the electron-proton system, the inflow ve- 
locity is O (0.1c a)- This agreement is unexpected 
since different mechanisms are clearly at play in 
the two cases — while there is some turbulence in 
the outflow jets of electron-proton reconnection, 
the amplitude is far smaller than in the pair sys- 
tem. However, since the growth of the Weibel in- 
stability appears to be crucial for maintaining fast 
reconnection in our pair plasma simulations, any 
suppression of the growth rate would serve to slow 
the process. One possible mechanism for doing so 
is the addition of a constant guide field component, 
B z . Such a field would tend to isotropize the tem- 
peratures in the current layer (T x w T y ) simply 
by keeping the particles on Larmor orbits and ei- 
ther substantially decrease or stabilize the growth 
of the instability. On the other hand, we also can- 
not rule out the possibility that another instabil- 
ity would develop when a guide field is present. It 
seems clear that any elongated current layer will 
ultimately develop strong pressure anisotropics, 
and in the presence of a guide field the parallel 
pressure p\\ will likely become much larger than 
that perpendicular to the magnetic field p±. For 
randomly oriented magnetic fields most reconnec- 
tion is guide field reconnection, so the anti-parallel 
geometry we consider here is in some sense a spe- 
cial case. But, if guide field reconnection is slow 
due to the suppression of the anisotropy instabil- 
ities, anti-parallel reconnection, even if rare, may 
dominate in pair plasmas. 

Even during anti-parallel reconnection the 
Weibel instability will be suppressed for large 
enough asymptotic plasma temperatures. As can 
be seen in Figure [H a large T can minimize the 
effects of the extra heating provided by the beam 
of outgoing plasma in the current layer. The 
critical temperature can be roughly estimated as 
T cr u ~ mv \l^ where va ~ ca is the drift velocity 
of the outflowing electrons. Thus for asymptotic 
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temperatures T > T crit « mc 2 A or f) > 1, the in- 
stability should be stabilized. Further simulations 
in this regime are being performed. 

After submitting this paper we became aware 
that Zenitani and Hesse also concluded that 
the Weibel instabilty plays a similar role dur- 
i ng reconnection in very r elativistic pair plasmas 

\Zenitani and Hesse\\200^j ) . 
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Fig. 1. — Reconnection overview. Panel (a): v ex 
over plotted with magnetic field lines. The solid 
and dashed portions of the lines indicate the direc- 
tion of the reconnecting magnetic field B x . Panel 
(b): B z at the same time. Panel (c): Cut of E y 
averaged horizontally over —50 < x < 50. Panel 
(d): Distribution of v y over the region —50 < x < 
50, — 2 < y < 2. Electrons are green, positrons 
blue. Panel (e): v z and v x distributions. v x is 
plotted in black and is essentially the same for 
both species. 
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Fig. 2. — Out-of-plane electron velocities for the 
simulations of Table HJ showing the evolution of 
the current layer with system size. Each panel 
shows about one-fourth of the total simulation do- 
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Fig. 3. — Reconnection rate versus time for the 
runs listed in Table [TJ The vertical lines indicate 
the times of the snapshots shown in Figure El Note 
the different scales on the horizontal axes. Differ- 
ences in the strength of the initial perturbation 
and the initial width of the current layer account 
for the differences in ramp-up times between the 
panels. 
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Fig. 4. — Panel (a): £^ at £ = 450, soon after 
the instability develops, for run c. Panel (b): Two 
components, T xx and T yy , of the positron temper- 
ature tensor along the line y = — 1. Panels (c)- 
(e): The velocity distributions taken between the 
dashed lines in panel (b) and — 1.5 < y < —0.5. 
Electrons are green, positrons blue. In panels (c) 
and (d) the two curves are nearly identical. 
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Fig. 5. — Panels (a)-(b): The structure of B z 
during the linear phase of the instability and af- 
ter saturation in an initially homogeneous plasma 
with T x /T y = 4. Panel (c): Theoretical (line) 
and measured (star) growth rate versus temper- 
ature anisotropy for the fastest growing insta- 
bility in a homogeneous plasma (finite k y with 
k x = 0). Panel (d): Growth rate in a Harris-sheet 
background. Stars indicate the measured values, 
the line the theoretical value in a homogeneous 
medium. 
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Fig. 6. — Out-of-plane electron velocity at St = 25 
intervals beginning at t = 600 for run d. 
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Table 1: Simulation parameters. 



Run Label 


Domain Size 


Gridpoints 


wo 


a 


100 x 50 


512 x 256 


2 


b 


200 x 100 


1024 x 512 


2 


c 


400 x 200 


2048 x 1024 


4 


d 


800 x 200 


4096 x 1024 


4 



Note. — wq is the initial half-width of the current sheet. 



Table 2: Simulation results. 



Run Label 




™out 


25 


2A 


*W 


^in,meas 


^in,calc 


a 


0.16 


0.27 


4.0 


35 


0.5 


0.13 


0.10 


b 


0.12 


0.32 


4.0 


80 


0.8 


0.15 


0.11 


c 


0.13 


0.33 


4.5 


120 


1.3 


0.16 


0.12 


d 


0.13 


0.30 


5.0 


135 


1.3 


0.13 


0.11 



Note. — 26 and 2 A denote the smaller and larger spatial extents, respectively, of the current sheet, fin meas and V[ n ca ic are 
the measured and calculated (from equation [5} velocities of plasma into the X-point. 
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